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(N '. Abstract 



We investigate the impact of the intrinsic heavy quark states on the predicted 



ON ■ values of CP asymmetries in the decays of heavy mesons. It is shown that 

the intrinsic charm contribution, although dynamically suppressed in QCD, 



is favored by the weak interaction, and therefore it can significantly dilute 
the predicted values of CP-violating asymmetries. This introduces additional 
non-perturbative uncertainty into the estimate of direct CP violating effects. 
We provide a phenomenological estimate of intrinsic charm content of r\ and 
rf mesons by expanding various amplitudes in terms of the heavy-light quark 
D . mixing angle and discuss theoretical uncertainties in the estimates of direct 



CP-violating asymmetries in B — > rf''K^*>. 
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I. INTRODUCTION 

Heavy quark decays offer a wide array of methods of testing the Standard Model and 
searches for the signatures of New Physics. In particular, they proved to be a powerful tool 
in the studies of the the weak mixing matrix and prepare a fertile ground for the exploration 
of the smallest CKM matrix elements, V u b and Vtd, and, most importantly, studies of CP 
violation. The observation of signatures of direct CP violation involves comparisons of the 
partial decay rates of charged mesons. As a result, the necessary requirement is a possibility 
for two distinct pathways (i.e. amplitudes with different weak and strong phases) to reach a 
given final state. It is realized, for instance, in the B decays to the final states not containing 
charmed quarks. In this example, CP violation can occur either from the interference of the 
tree-level and penguin amplitudes, or from the interference of the penguin amplitudes with 
different quark flavors in the loop. The strong phases are generated by allowing internal 
quarks in the penguin loop to go on their mass shells by virtue of the so-called Bander- 
Silverman-Soni (BSS) mechanism ||1|. It is expected that the resulting strong phases for 
different amplitudes are different, in one case because tree level graphs do not produce 
perturbative strong phases, and in the other, because of the different mass thresholds for 
different quark species in the penguin loop. 

Let us consider a process governed by the quark b — > s transition to a given final state /. 
In general, there are two distinct ways to reach /: by the tree- level b — > uus amplitude or by 
a penguin b — > sqq transition. It is clear that the tree-level amplitude is proportional to the 
CKM matrix elements V* b V us and thus Cabibbo-suppressed as it scales as A 4 in Wolfenstein 
parameterization. On the other hand, the leading penguin effects, although suppressed by 
loop factors, scale as V t lV ts , or A 2 . This makes the tree amplitude comparable in strength 
with the one-loop penguin diagrams, thus enhancing the interference term and allowing for 
sizable CP-violating asymmetry 

A% = l^'ll*" - (1) 

These asymmetries have been calculated for a number of final states /, such as B — > 
KK, Kir, etc. 

It is plausible to assume that while BSS mechanism does represent a way to produce a 
non-zero CP violating asymmetry, the soft, non-perturbative effects might produce larger 
FSI phases in the exclusive transitions, thus introducing a non-perturbative uncertainty to 
the calculation of the asymmetry of Eq. ([!]). For example, the rescattering of physical 
hadrons produced in the reaction provides an additional source for the FSI phase @-§|. 

Here we shall argue that soft FSI contributions do not exhaust list of possible non- 
perturbative uncertainties of Aq p . The Fock state expansion could include a non- vanishing 
contribution of the heavy quark (e.g. charm) states to light meson's wave function. Although 
higher Fock state contributions are dynamically suppressed in QCD, weak transitions of the 
b quark to the heavy quark states are Cabibbo-favored. Therefore, weak interaction selects 
intrinsic charm states of the light mesons making their contributions competitive with the 
direct b —>■ u transitions to the leading Fock states of the light mesons. We admit that the 
intrinsic charm content of light mesons is considerably difficult to estimate. However, at 



least for a particular class of light vector and pseudoscalar mesons this contribution can be 
phenomenologically accounted for by allowing "mixing" of the heavy mesons with hidden 
charm with the light mesons bearing the same quantum numbers. For instance, there is a 
non-vanishing probability for the mixing of J/ip and <fi, t] c and 77', etc. In what follows we 
consider the possible effects of intrinsic heavy quark states of light mesons on CP violating 
asymmetries. 

The paper is organized as follows. In Section II we consider an upper bound on the 
value of the heavy-light quark mixing angle. In Section III we discuss how intrinsic charm 
content of the light mesons affects direct CP violating asymmetries concentrating on the 
phenomenologically interesting transitions B —* rf^K^*\ We summarize our results in 
Section IV. 



II. HEAVY-LIGHT QUARK MIXING 

There are many possible approaches in QCD that account for the intrinsic heavy quark 
states in light quark systems. It is therefore reasonable to employ a phenomenological 
description of the mixing, extracting the values of the mixing angles from the experiment. 
These values can be later compared to the results obtained using various models. Let us 
parameterize a mixing of the heavy and light pseudoscalar mesons in terms of the following 
matrix 

ap — bp cp — dp — sin ap cos(0p — 9p) \ { Vo \ 

—cp — dp ap + bp —sinapsm((j)p — 6p) | r] 8 , (2) 

sin ap cos 4>p sin ap sin <j) P cos ap J \ r] c0 J 

where ap = cos 2 (ap/2) cos9 P , bp = sin 2 («p/2) cos(20p — Op), cp = cos 2 (ap/2) sinO P , and 
dp = sin 2 («p/2) sin(20p — Op). Here 770, 77s, and r/ c0 represent the flavor SU(3) singlet, octet 
and pure cc states respectively, and 77', r\, and r\ c are physical states: 

1 770) = —j=\uu + dd + ss), 
V3 

\r] 8 ) = —=\uu + dd — 2ss), (3) 

V6 

|?7co) = |cc). 

The mixing angles ap, Op, and <pp have a simple physical meaning 0: ap represents an 
"admixture" of the heavy quarks to the light ones, Op represents a mixing of the light quarks 
among themselves, and (f)p gives a light quark 577(3) admixture to a heavy QQ state. It 
is clear that one recovers a standard 77 — 77' mixing matrix as ap — > and heavy quarks 
decouple from the light ones. Sometimes, it is more convenient to work in the quark basisf] 




*Note that we are not using any specific quark model. 



\f]') — —7=X v t\uu + dd) + Y v t\ss) + Z v i\cc), 
V2 

\r]) = —=X v \uu + dd) + Y v \ss) + Z v \cc), (4) 

V2 

\r] c ) = -=X Vc \uu + dd) + Y Vc \ss) + Z Vc \cc), 

where we have generalized the construction of to include intrinsic charm states. Using 
Eq. dD it is easy to show that 

K + K + Z l = l ( 5 ) 

for each r\i = {r] , rj' , rj c } . This equation might be violated by the presence of some other 
pseudoscalars that might mix with r\i as well. This would result in the limitation of the 
described method to the prediction of the upper bound on the value of r] c contribution to 
the mixing angle ap. 

The construction of Eq. (||) corresponds to the introduction of the additional, charmed 
singlet current in addition to the "standard" SU(3) singlet and octet currents 

A% = -j={uYl 5 u + d^^d - 2sYl5s), 

A% = -j= {uYl5U + dYlsd + sYl5s), (6) 

These induce the following matrix elements 



W\uj^5u\0) = -i\h ( a P ~ b p) F o + ~7k( c p ~ d p) F s 




L V^ J 

(»/|s7mT6«|0> - -i\\- [(a P - b P )F - V2(c P - d P )F 8 ] p M , (7) 

(V|c7m75c|0) = iv / 2sina P cos(0p - 0p)F VcO p tl . 

Similar matrix elements exist for the r\. Here F 0t8 are the singlet and octet decay constants 
and F Vc0 is r/ c o decay constant. In the limit of SU(3) symmetry F 8 = F n . The 5f/(3)-violating 
corrections have been calculated in || in the framework of Chiral Perturbation Theory and 
found to modify this relation by approximately 25%, i.e. Fg/F n = l/x$ = 1.25. We shall 
use this value in the following analysis. The value of Xo, on the other hand, is not fixed by 
the SU(3) symmetry arguments, so we shall keep it as a free parameter fixing it later by 
fitting to the experimental data (In the limit of the "nonet" symmetry Fo/F n = 1/xq = 1). 
The parameters of the mixing matrix can be obtained phenomenologically, as they con- 
tribute to the decays of charmonia to the light 77, 77' and 7r° mesons, and to radiative decays 
of the light mesons. Here our essential assumption is that the mixing angles ap and <fip are 
sufficiently small. This assumption, however, is rather loose, and is valid even for a relatively 
large charm content of 7/ |9|{ll|) but it allows us to carry out a perturbative expansion in 
these angles. In the following analysis we will only keep terms linear in ap and (pp. 



The bulk of information about the relevant mixing angles comes from the radiative decays 
of r), 7]', and t] c mesons. Normalizing the decay widths to the width of tt° — > 77 and using 
Eq. (El) we obtain 



m, 



/V 



r(V -» 77) 

\m v/ J r(7r°^77) 

F n F^ c P — dp 

— [a P - b P ] + — - 

-TO -T8 



l 3CF w F r]c0 



in 



sinapcos(0p — 0j 



<lc 



Pv 



m Tl 



m r 



Y( v -» 77) 

r(7r° — >• 77) 



1 CFirF Vc0 



^[a P + bp ]-^V8{c P + d P ] ^ 



sin ap sin(0p — 9j 



PVc 



m T , 



m r 



F(r]c -> 77) 

r(7r° — > 77) 



F n . F n sin a P sin P \ZCF^F ncD 

— sm a p cos <f>p + — = \- \ — — cos ap 







m. 



<lc 



where ap, b P , and cp were defined in Eq.(0) and 

Expanding in terms of ap and (f)p and keeping only the linear part we arrive at 

2 
2^8 . „ „ / /3CF 7r F„^ n \ 

A?' 



(9) 



—j= sm c^p + cos 0p \Xq — \ — 7, dp 



in. 



P>i 



Pvc 



x 8 cos9 P - sin9p V8x - \ n 2 — a P 



(10) 



x a P + 



13 Cr n r Vc0 



711. 



>lc 



As it is seen from Eq.(|10l), the dependence on <pp drops out at this order, so we set <pp = 
in what follows. In time, when the accuracy of experimental measurements improves, the 
second order in the "angle expansion" should constrain the value of <pp as well. 

It is clear from Eq.(|T0"D that three equations do not constrain four parameters: 6p, ap, x , 
and F Va0 . There are several ways to proceed at this point. For instance, one can fixio = l 
by assuming a nonet symmetry and fit the rest of the parameters from the Eq. fllOD . Here we 
shall take a different approach. The "angle expansion" can be carried out for other processes 
involving 77 and 7]', such as J/ip — > 77^7 and J/ip — > 7/ c 7, or rj — > fry. In the limit of SU(3) 
invariance, only the singlet state 770 or 7] c0 can be coupled to J/if). The amplitudes of the 



radiative decay of a vector charmonium state into a final state containing the SU(3) singlet 
light quark state \rj ) or |t7 c0 ) can be written as 



A(J/ip -> r]o7) = Ae^c^e^cp^h . 



Ae^ape^e^p^k^ 
A{J/^ - rial) = Be^qfyffi 

Performing the "angle expansion" we find that [] 

r(j/v> - ^7) 1 



(ii) 






7/7) 
^7) 



cos 2 6 1 



P \_P n > _ 



Y{J/i, 
T{J/*l> 



777) 
7/7) 



sin 2 #. 



r(J/^ -> 777) tan 2 # P 



.#•7. 

A/ 
Pt7 



l3r 



l + q P (A/fl) 
ap - (A/5) 

1 + a P (A/B) 
a P - (A/B) 



n 2 



n 2 



(12) 



As seen from Eq. ([12]) , the ratio of radiative decay widths of charmonia into 77' and 77 is 
independent of the mixing angle atp at this order and can be used to extract the value of 
dp. Alternatively, the ratio 



T(r/ / -» fry) 
3T(uj — > 7T7) 






A^, 



1 r /- " - 

- v2cos^p + sin^p 



(13) 



does not depend on the heavy-light mixing angle and can be used to extract 9p. Similar 
analysis is possible for p — *■ 777 etc. Extracting #p from either decay mode and feeding it 
into Eq. fllPf) leads to the constraints on all of the mixing parameters. This calculation finds 



0, 



-20°, £0 ~ 0.92 (which corresponds to Fq/F^ fh 1.05 in full accord with 



The 



situation is more complicated with respect to the heavy-light mixing angles. Uncertainties in 
the extractions of the light quark mixing angle Op and decay rates complicate the extraction 
of the value of the mixing angle a. p. In fact, all of the experimental results can be successfully 
fit assuming a zero value for the mixing angle a. p. In order to overcome this problem, 
additional constraints on the parameters in Eq.(^) must be imposed. For instance, the 
value of F Vc0 can be extracted from other decay modes which are less sensitive to the heavy- 
light mixing, e.g. B —>■ rj c X s , where the cc component is enhanced by the CKM matrix 
element V^. Using this method, F Vc0 can be fixed to F Vc0 ~ 0.29 GeV. This and the third 
line of Eq.QTOD implies that | sinap| < 0.03 (1.7°), <pp = 0. It is interesting to compare this 
bound to what already exists in the literature. It appears that all the models of intrinsic 
charm based on heavy-light meson mixing [|T1|,|T^] (see also [I3~|) satisfy this bound. On the 
other hand, there exists a class of OPE-based calculations P,|10|| that clearly violates this 



^This result agrees with the result of [fj [| in the limit A/B <C ap. Please note that it is not possible 
to extract ap from this decay mode alone without invoking additional dynamical arguments about 
the size of A/B [cf. @]. 



bound, predicting ap ~ 7° if the expansion is terminated at the level l/m 2 c . As we shall see 
later, these values of ap significantly dilute the direct CP asymmetries in B-decaysQ. 

A similar construction is available for the vector mesons, where we use the same notations 
with the obvious replacement of the subscript P by V 



(14) 
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( ay — by 


Cy — dy 


— sin ay cos((py 


-0v)\ 


( Vq 


u, = 


—Cy — dy 


ay + by 


— sin ay sin(0y 


-Ov) 


L0$ 


J It) 


\ Sin ay COS (j)y 


sin ay sin (f) V 


COS tty 
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Uo 



Numerically, sin ay cos( 



7 v 



1.2 x 10" 



III. AMPLITUDES AND CP- VIOLATING ASYMMETRIES 

In fi-decays, the tree level amplitude is suppressed by a small V u b, which makes it 
comparable to the one-loop penguin amplitude. On the other hand, the transition B — > 
[QQ]s — > Ms with Q being a heavy (charm) quark and M a light final state meson is 
not CKM suppressed. Therefore the latter mechanism becomes competitive with the direct 
b —>■ u transition to the light quarks constituting the light mesons. Since b — > ccs amplitude 
does not contain a CP violating weak phase, it may significantly reduce the predicted value 
of CP asymmetry. In what follows we consider two cases for M being a pseudoscalar and a 
vector meson. 

The generic amplitude for the decays of a B meson to the charmless final state / can be 
written as 



A B ^ f = £ U A T + £ C A M + ]T i t A P = £ u \A T + Af\ + £ 



i=u,c,t 



A c ' + A 



M 



(15) 



Here At is a tree-level amplitude, Am is a "mixing" amplitude and Ap is a penguin ampli- 
tude. As usual, unitarity of the CKM matrix was used to write £t = — £ c — £„, Ap = Ap — A P , 
A P l = A U P - A P , and & = V* b V is . 

In order to form a CP-violating asymmetry, we also must consider the corresponding 
amplitude for the decay of the B, 



At 



B^f = C A T + A? + C 



A% + A 



M 



(16) 



Using the fact that At = Ap and Am = Am, the asymmetry of Eq. ([]]) can be formed as 



^The OPE-based calculations successfully explain the unexpectedly large branching ratio for 
B — > ij'K. They however have certain phenomenological difficulties in the case of the inclusive 
rf production in B-decays, as well as describing ratios of branching fractions of B meson decaying 



to PP and PV final states [11,26,27| . The fact that the mass of the charmed quark is not suffi- 
ciently large for the fast convergence of the OPE might explain the variance in the results of these 



calculations pA \. 



A, 



r 



B->f 



B-+f, 



A f = X f Im£*£ c Im 



A T + A 



ut' ! 



At + A M 



(17) 



where A/ = J (I — (x/ x + Xf 2 ) 2 )(l — (xf x — s/ 2 ) 2 )/(47rms) is a phase space factor, divided 
by the sum of the corresponding decay rates. As seen from Eq . (|17D , direct CP violation in 
the B — > T]iK mode arise not only from the interference of the penguin diagrams with the 
different internal quark flavors, but also, due to the complicated quark content of the rji, 
from the interference of the Cabibbo-suppressed tree-level amplitudes with the penguin am- 
plitudes. As we shall see later, intrinsic charm contribution affects CP violating asymmetries 
calculated using perturbative BSS phases. It is therefore instructive to study the dependence 
of the asymmetry on the parameter q = q 2 /m 2 B which parameterizes the momentum flowing 
through the penguin vertex. The final answer, of course, is independent of q and should be 
obtained by either integrating the asymmetry with respect to q smeared by some function 
defined by the momentum distribution of quarks in the final state mesons, or by fixing q us- 
ing quark model arguments. We shall use the second method for our predictions, presenting 
the graphs asymmetry vs. q to show the threshold structure of CP violating asymmetries. 
We comment on the effects of soft FSI phases in the Conclusion. 

In the following discussion we shall first use the effective Hamiltonian calculated at 
leading order in QCD, i.e. with no QCD corrections associated with the penguin part. 
Next, the full Next-to-Leading order effective Hamiltonian is employed. 

It is well-known that the calculation of the two-body non-leptonic decays of heavy mesons 
cannot be performed without invoking a particular model. This model dependence is par- 
tially cancelled in the asymmetry Eq.(|l|). In our calculation we choose factorization approx- 
imation |I7-2I| and the Bauer, Stech and Wirbel (BSW) model |15[ to estimate relevant 
form-factors. 



A. Leading order calculation 



The "no QCD Hamiltonian" reads 



njLO 
H eff 



^(4<**-5F 



i=l 



£ tfi 



i=u,c,t 



° 3 +o ° 5 +o 



(18) 



01 = Sa^LQpQ/ffLba , 0\ = s^LQQ-fLb , 
3 (5) = 87/i* Lb lLin' iI '(R)Q > °4(6) = s a -i^Lbp^2q^L{R)q a 



Here, Q = {u,c}, q = {u,d,s}, and Fi are the Inami-Lim functions for the flavor i. A 
similar construction is available for the effective b — ► d transitions (although the effects of 
intrinsic charm states are largely suppressed in these modes since b — > uud decays are not 
CKM suppressed compared to b — ► ccd). In what follows we drop the contributions from 
the electroweak penguins and dipole operators for the sake of simplicity. 
The decay amplitude A (,) K = —i(r]^K\TC e ff\B) can be written as 



A 



B->r](')K 



Zu 






l\rp 



(>) K 



T 
—GpTT? B 

a 



ut-^-p 



+ £ c 



F ct Af K + Af K 






Af« = G F ml^ I 1 



a2^f+ K W )^(^) + aiF K f?l (p)Lr,(jii 
1 \ 



(19) 



m 



a P , 



A V M K = G F m 2 B sma P cos(4> P - 9 P )F V J^ (m 2 {/) )L k (fii), 
where d\ = C 2 + X^i, a 2 = C\ + X&2, Hi = "^f/^Bj and the following notations are used: 
a P = F K ff\0)L v (jii) + F'M?(< lt) )LM 



+ 2F K ff(0)^^M v (x 2 ) + F^f^{^)-^-m b M h (jM). 

rn rn, ' ' nm 



(20) 
m s mb ' " " ''" ' ''" m s 

Our choice of the form factors and kinematic parameters is explained in the Appendix. 
In the standard factorization approach \ = 1/N C and all the octet-octet non-factorizable 
corrections are neglected. These corrections can be accounted for phenomenologically, by 
treating \ as a free parameter and fitting it to the available data assuming universality of 



these corrections for different final states p2[ . In this approach, a 2 — 0.25 p3| , |24| . 

In order to maintain gauge invariance and unitarity, the calculation must be performed 
to order a 2 in the PQCD expansion |[25|| . The CP violating asymmetry reads 



\(')K = \(0K Im CJ,c 



Af K ) 2 [ 



ImF*ReF ct + ReF, 



Lit 



ImFat - ^ReF ct 
o 



+ Af K Af K ImF* ut + Af K 4' )K 



tip 
ImFct - —ReF ct 
o 



(21) 



Dividing A v 'k by the sum of the decay rates, the asymmetry Eq. (|J) can be calculated. As 
it is seen from Fig. (|l|), intrinsic charm reduces CP asymmetry by approximately 30 — 50%, 
if our estimate of a>p is used thus complicating the extraction of CP violating parameters of 
the CKM matrix from the decay modes of this type. Please note, that if by some reason, 
the intrinsic charm content of the rf is increased P,|TU||, this mode becomes practically 
useless for the observation of the direct CP violating effects. It is clear from the Fig.(^) 
that CP violating asymmetry is significantly diluted even for otp ~ 7° which corresponds 
to the lower bound of the prediction in OPE-based calculations. This however, simplifies 
the time-dependent analysis of the decay B® — > f]'K s as it suppresses direct CP violating 
amplitude. The behavior of CP-violating asymmetries for the case of r\K final states is 
similar to the described above. We therefore refrain from displaying the shape of A^ P (q) 
function, but rather show the numerical value of Aq p in Table 1 with q fixed by the quark 
model arguments. The analysis of the rp''V final state is completely similar to the calculation 
described above. The results are presented in Table 1 as well. 



B. Next-to-Leading order Calculation 

The NLO QCD effective Hamiltonian reads 



9 



n/NLO _ ^^F \ sr^ C 



6 



^CA g + E^o 



Li=l fc=3 



(22) 



0? = s^LQpQ^Lb,, , Of = s^LQQ^Lb , 
03(5) = s^^Lb^q^L(R)q , 4(6) = s a ^^Lbp^q^L(R)q a . 

In our calculation we used renormalization scheme independent effective Wilson's coefficitens 
defined as 

C 2i+i Q ( m b) = C 2i+ i(m b ) + g * f G(m g , g 2 , w 6 ) - yj C 2 , 

<#' 9 K) = C 2i (m b ) - ^p^ (G(m q ,q 2 ,m b ) - H) C 2 (23) 

with i — 1,2 and Wilson's coefficients C& at /i = rrib with a s (m^) = 0.118 are given by 

Ci(ra b ) = -0.313, C 2 (m b ) = 1.150, 
(73(m 6 ) = 0.017, C 4 (m b ) = -0.037, 
C 5 (m b ) = 0.010, (7 6 (m 6 ) = -0.046. (24) 

The function G{m q ,q 2 , //) is given by 

2 /■ rrig — x(l — x)q 2 

G(m q , q , fj,) — —4 / dx x(l — x) log — (25) 

J fi 2 

where q is a momentum of the gluon in the penguin diagram. It is clear that the strong 
phase is generated every time q 2 > 4m 2 for each quark species in the loop. 

The factorization calculation is completely similar to the one performed in the previous 
section. The only difference arises from the fact that the effective constants a^ = C^-i + xC^ 
are not vanishing and have to be taken into account along with the corresponding form- 
factors. Those are the combinations of the form-factors and the decay constants defined 
in the previous section. The results are presented in Table 1. In our calculation we fixed 
the value of a 2 — 0.25 fixed by the experimental data while dropping the nonfactorizable 
contributions to the matrix elements of penguin operators which are difficult to estimate 
reliably. The asymmetry is seen not to change significantly in going to NLO approximation 
with the intrinsic charm contribution still diluting it at the level of 30-50 % for the estimated 
value of ap ~ 0.03. The effect becomes stronger for higher values of mixing angles. As it 
is seen, even the perturbative result is very uncertain. The major source of uncertainty 
is by far dominated by the value of the charmed quark mass that affects the position of 
the charmed quark threshold. Because of the 577(3) symmetry relations, the values of the 
hadronic form-factors, although important for the decay width predictions, are not seen 
to significantly affect the predicted values of CP asymmetries. An additional source of 
uncertainty is the value of non-factorizable corrections usually summarized in the effective 
parameter \- I n the case of the asymmetries induced by the penguin-tree interference, it 
can change the balance of these contributions, thus shifting the value of Aq P . It is usually 
assumed that the non-factorizable corrections can be taken into account by replacing the 

10 



factors of 1/N C in Oj throughout the calculation with the single parameter \. This fact also 
induces the uncertainty into the estimate of both decay width and CP- violating asymmetry.| 
It is interesting to note that in the limit x ~ > the contribution from the intrinsic charm 
amplitude becomes negative actually reducing the predicted values of branching fractions: 
for instance, Br(B~ — > T]'K~) = 2.8 x 1CT 5 for moderate values of the form-factors and 
ap ~ 0.3 drops to Br(B~ — > rj'K~) = 1.6 x 10~ 5 if x — 0. Of course, higher values for the 
branching fraction are still possible considering large uncertainties in the values of hadronic 



form-factors 26 



Similar calculation can be performed for the decays involving vector particles in the 
final state, such as B —>■ <pK. The calculation is simplified considerably since the decay is 
dominated by a single penguin amplitude. The CP asymmetry for the leading order case is 



A#r = X^k PmCuic 



Af 2 {ImF: t ReF ct + ReF^ImF*} + A^ A P K JmF* 



(26) 



Here, the following notations are used 

A P K = -G F m^{p B +p K ) ■ e ff ' 

a< m = -Gfui^F^pb+pk) ■ e f+(ml)a 2 sin a v cos(0 y - 9 V ) (27) 



57T 



As usual, one expects partial cancelations among the first and the second term in Eq.(p6|) 
because of the GIM mechanism. Thus, the intrinsic charm amplitude is potentially important 
as it does not suffer from this cancellation. Fortunately, the small value of the heavy-light 
mixing angle for the vector mesons makes the intrinsic charm contribution extremely small. 
The correction to the asymmetry Aq P ~ 0.3 % is less than 1 % for the same choice of CKM 
parameters as before. We therefore see no point in performing the NLO studies of this 
class of decay modes. Clearly, final states containing vector particles are much less affected 
by higher Fock state "pollution". The resulting CP violating asymmetries, however, are 
significantly smaller. 

IV. CONCLUSIONS AND OUTLOOK 

We have investigated the impact of the intrinsic heavy quark states on the CP violating 
asymmetries in B-decays. It arises due to the fact that the intrinsic charm quark states, 
although suppressed by QCD dynamics of the process, are less affected by weak Cabibbo 
suppression. The most dramatic effect is seen to occur in the case of the rfV final states. 
Unfortunately, the impact of the intrinsic charm states on the CP asymmetry is very difficult 



§The contribution of the octet matrix element is usually associated with the quark final state 
rescatterings into the colorless hadron states. It is clear that the interactions of the rescattered 
and spectator quarks might introduce different x' s f° r different amplitudes. For instance, there 
would be different contributions from b — ► uus and b — ► dds octet amplitudes to B~ — ► rft'K^*'. 
The effect is similar to the effect of Pauli interference in nuclear physics. 
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to test experimentally since the direct CP asymmetry explicitly depends on the strong phase 
of ^-matrix element relating two states produced by the weak interaction. This quantity 
is notoriously difficult to estimate theoretically since it comprises not only perturbative 
BSS phases but also soft, nonperturbative phases generated by the rescattering of physical 
hadrons 0. Soft FSI contributions, although important, are not seen to change the shape 
of the graphs asymmetry vs. parameter q but rather move them up or down, whereas the 
described mechanism certainly affects the shape. Of course, in the final result the function 
described by Afj r P {q) has to be integrated over q with a suitable choice of "smearing function" 
representing momentum distribution of quarks inside the final hadrons, or q has to be fixed 
by the quark model arguments. In either case, the intrinsic charm contribution is seen 
to reduce the value for CP violating asymmetry sizably. This effect is not universal, but 
rather specific to the final states containing rf 1 ' mesons, so it cannot be attributed to the 
long-distance part of the penguin diagram. 

The mechanism described in this paper, along with the uncertainty in the determination 
of FSI phases complicates the use of the direct CP violating asymmetries as a so-called 
"consistency check" in the determination of the angles of the CKM triangle. For instance, the 
numerical values of the angle 7 determined from various decays modes have to be consistent 
with each other unless there exists a New Physics contribution that violates this requirement. 
Thus, the possible inconsistency should manifest a Non-Standard Model mechanism affecting 
the decay processes. As one can see from the discussion above, there exist possible sources 
of violation of this consistency check within the Standard Model. 

Finally, we would like to note that this mechanism does not markedly affect the decay 
modes where electroweak (EW) penguins are manifest (e.g. B — > rj'ir) since there the tree- 
level and mixing amplitudes contribute at the same order in Wolfenstein parameter A with 
mixing amplitudes additionally suppressed by the small values of heavy-light mixing angles. 

ACKNOWLEDGMENTS 

I would like to thank Gene Golowich, Barry Holstein, Nathan Isgur, Alex Kagan, Amarjit 
Soni, Daniel Wyler and especially John Donoghue and Adam Falk for useful discussions and 
insightful comments. This work was supported by the National Science Foundation under 
Grants Nos. PHY-9404057 and PHY-9457916, and by the Department of Energy Grant No. 
DE-FG02-94ER40869. 



APPENDIX A: HADRONIC MATRIX ELEMENTS AND KINEMATICAL 

FACTORS 

We use the following definitions for the pseudoscalar hadronic form-factors 



(M\qj,b\B) = f^(q 2 )(p B +PM), + /5(g 2 )(p B - Vu, 



jU) 



(M\qu,j 5 q 2 \0) = -iV2F^p Mll . (Al) 

We take Fk = 0.12 GeV. Usual relations among decay constants (and transition form- 
factors) of different mesons in the pseudoscalar octet imposed by S'6 r (3)-symmetry are mod- 
ified in the presence of r\ — rj' mixing and intrinsic charm components of f] and 7/, e.g. 
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TPUU 



-y=Ua P - b P )F + —={c P - d P )F 8 



F°? = -^ {(a P - b P )F - V2(c P - d P )F 8 } , 



(A2) 



where F ~ 1.05F n ,F 8 ~ 1.25F n as explained in the text, and F w = 0.093 GeV. Also, 
/° = 0.33. Similar expressions exist for the ^-meson as well. Vector and axial vector 
form-factors are defined as 



2\„* 



(1^1^(1+75)615) = ig{q 2 )e, vaP e* u {p B +p K ) a (pB ~Pk? + (m% - m 2 K )h{q')e 
+ (e* • q) [f2(q 2 )(PB+PK)t, + fs(q 2 )(PB - Px) M 
(M*|gi7 M g 2 |0) = -iV2F M *m M *e* 



(A3) 



The relationship of the form-factors used in this paper to those of [[0| is 



S(<7 2 



iV(q 2 



m B + m/f 



, fi(q 2 



iA t {q 2 



nip — rriK* 



hW 



iA 2 {q 2 



m B + m K * 



h{q 2 ) = -o [2m K *A - (m B + m K »)A 1 + (m B - m K *)A 2 ] 



(A4) 



where we take A\ = A 2 = A% = 0.33 GeV for the sake of simplicity. L(/j,i) and M are the 
kinematical parameters: 

■L>k\Jii) - 1 - U-K + , K / 2)^1 

T ( \-1 fV ~ {q2) 

- L r){(J'i) — J- — Hi) ~r Ttfi 2 \^K, 



M k (ni 



(3 - y + (1 - 3y)fi K ~ (1 - y)lb, 



f K (q 2 ) 



■ " K r 2 a i - y + (i - v)iiK + (i + yK) 



(A5) 



and y ~ 1/2 + 1 is related to the momentum distribution of quarks inside of the mesons. 
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TABLES 



in, 



TABLE I 

= 1.6-5-1.3 GrV 



CP- Violating Asymmetries Afip for q 



q 2 jm\ ~ 0.5, p = 0.05, rj = 0.36, and 



Mode 


4g£(ap = 0.0), % 


Af P {a P = 0.03), % 


A d j r P (a P = 0.12), % 


B~ ->r/K~, LO 


1.9-r3.5 


1.6-5-2.7 


l.O-i-1.5 


B- ->r]'K-, NLO 


2.5-5-4.4 


1.7-5-2.9 


0.7-5-1.2 


B- -+T)K~, LO 


3.1-7-7.3 


2.7-5-6.1 


1.9-5-3.9 


B~ -+r)K~, NLO 


3.0-7-7.1 


2.2-5-5.4 


0.8-5-2.6 


B- ->t/K*~, LO 


2.9-5-6.7 


1.9-5-4.1 


0.9-5-1.6 


B~ ->t/K*~, NLO 


7.6-5-16.5 


3.3-5-7.5 


0.3-5-1.3 


B- ^rjK*~, LO 


6.7-5-19.8 


4.5-5-12.5 


1.9-5-4.7 


5- ^rjK*~, NLO 


15.4-5-38.3 


6.3-5-18.3 


1.4-5-1.9 
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FIGURES 




FIG. 1. CP asymmetry calculated at the leading order in QCD as a function of the parameter 
"'" " ,J '"'■ n — ^ T]'K, without intrinsic charm (grey) and with intrinsic charm ap = 0.03 (black). 



1 /m 2 B for B 
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FIG. 2. Same as Fig.(l) but for a P ~ 0.12 . 
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